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Abstract 

> 
On 

We calculate limits on the trilinear soft-breaking parameter, At, in the Min- 

\Q 

imal Supersymmetric Standard model by requiring the absence of nonzero 

in 

Q\ ' top squark vacuum expectation values. Assuming a low tan f3, which implies 

Q_i' a large top Yukawa coupling, we also calculate one-loop corrections to the 

effective potential. The resulting numerical calculations of the charge/color 
breaking limits are presented as best-fit surfaces. We compare these results 
^ ' with the analytical limit, A\ < 3(to2 + ml + mf), and find that although this 



q > ■■•'t- 
is a good estimate of the charge/color breaking bounds for a simplified model 

of the top sector, stricter bounds are found by a numerical minimization of 

the Minimal Supersymmetric Standard Model potential. 

I. INTRODUCTION 

The Minimal Supersymmetric Standard Model (MSSM) may be described by a La- 
jrangian containing interactions consistent with invariance under the gauge group SU(3) x 
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SU(2) x U{1)y and global supersymmetry plus a Lagrangian containing a restricted set 
of soft supersymmetry breaking terms These terms break supersymmetry while main- 
taining a useful property of a supersymmetric theory, namely the cancellation of quadratic 
divergences 0. The absence of these divergences is necessary in order to define the renor- 
malized mass of a fundamental scalar, such as the Higgs boson, without a fine-tuning of the 
cancellation between the bare mass and the scalar self-energy ||. 

The presence of fundamental scalar fields in the MSSM, besides the Higgs bosons, leads 
to the possibility that these fields may acquire non-zero vacuum expectation values (vevs). 
Since this would violate the conservation of color and/or electric charge symmetry, this 
leads to forbidden regions of the parameter space of the theory. We will calculate numerical 
estimates of the boundary of the allowed region of soft-breaking parameters using both the 
tree-level potential and the one-loop effective potential. 

Many studies of the MSSM mass spectrum neglect these charge/color breaking, or CCB, 
bounds in their analyses. Previously, CCB bounds were obtained for various supersym- 
metric models, however no systematic numerical study of CCB constraints for a realistic 
approximation to the MSSM using the one-loop effective potential has been done HJ5J. 

One may assume that there are relations among the soft breaking terms, such as in the 
minimal supergravity model in which all scalar masses and scalar trilinear couplings are 
the same at the unification scale, of order 10 16 GeV 0. However we will find constraints 
on the soft-breaking parameters at a low-energy scale, Qq, with Q < 1 TeV. This is an 
indeterminate upper limit on particle masses if the MSSM is to explain the gauge hierarchy 
problem. We will not make any assumptions about the theory near the GUT scale nor the 
particle spectrum above Qq. 

We will use an approximation to the MSSM that includes only the top flavor supermul- 
tiplets. This follows from evidence that the top quark mass m t « 176 GeV 0. We use 
the conventional definition tan/3 = i^/Vl, with v\, v% the vevs for the Higgs scalar fields, 
Hi and H 2 , respectively. Assuming a small value for tan/3, near 1.0, gives the top quark 
Yukawa coupling, h t = 1.0. The contributions from the bottom supermultiplets may then 
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be ignored. 

There are various reasons to choose these particular values of tan/5 and to consider 
only the top squarks as acquiring a non-zero vev. First of all, there is an infrared quasi- 
fixed point in the renormalization group equation for h t (Q) which corresponds to a value 
hf p (m t ) ~ 1.25 || . The mass relation 

h FP (m t ) 

m t (m t ) = p vsmp (1) 

r(2) 



gives tan/3 pa 1.2 if one uses the relation between the top quark mass defined by the pole in 
its propagator and its running mass to first-order in the QCD gauge coupling ||. 

4 

m po i e = mt(m t ) 1 + — a 3 (m t ) + 0(af) . (2) 

Therefore a value of tan/3 at Qo in the range 1.1 < tan/? < 2.0 results from a large range of 
h t values at the GUT scale. Although tan (3 is not required to be in this range, it indicates 
that this is a natural choice. One motivation for considering only the top sector comes 
from assuming common soft-breaking parameters at the GUT scale. A large value of h t 
causes the third generation parameters to undergo the largest change as they are evolved 
from Qgut down to Q . For this same reason, h t also gives the largest contribution to the 
radiative gauge symmetry breaking of SU(2) x U(1)y [pXi| , [TT| . Therefore, if one assumes 



that the minimum of the effective potential at energy scales Q 3> Qo gives zero vevs for 
the scalar fields, such as in the case of universality at the GUT scale, as one evolves to 
Qo the third-generation parameters undergo the largest change and the CCB constraints 
from third generation scalar fields will be the most restrictive. Finally, as discussed in 



Ref. ||12|| , the potential barrier height for tunneling from the the symmetric vacuum at a 
high temperature (T > 100 GeV), early in the expansion of the universe, to a lower CCB 
minimum is proportional to 1/ h? min where h min is the smallest of the Yukawa couplings for 
the slepton and squark fields that have non-zero vevs at the CCB minimum. This implies 
that one should consider CCB vacua in which only the Higgs fields and the top squarks have 
non-zero vevs in order for the tunneling from the symmetric to the CCB vacuum to have 
occurred in a time less than the present age of the universe (~ 10 10 years). 



II. APPROXIMATION TO THE MSSM 



We use a consistent approximation to the MSSM with tan/3 = 1.5 as a small value 
near the fixed point value and interactions with the bottom quark superfields are ignored. 
We use all MSSM interactions between the following fields H±, H 2 , Hi, H 2 , q, t c , q, t c , 
A^, and A. Hi and H 2 are respectively the hypercharge —1 and +1 Higgs boson doublets. 
The corresponding field variables with a tilde are the Higgsino doublets, q and t c are the 
left-handed component of the top quark and the right-handed component of the charge 
conjugate top quark field respectively. Again, the corresponding field variables with tildes 
are the top squarks. is the gluon field and A is the gluino field. Notice that the field 
content in this approximation is supersymmetric. This arises from including all interactions 
with the top quark supermultiplet involving the parameters h t and fi. The potential in this 
approximation as well as the definitions of the parameters appearing in it are shown in the 
appendix. 

We use the values of the gauge couplings at the weak scale Mz- g\ = 0.358, g 2 = 0.652, 
and gs = 1.213 for all calculations [jH5| |. Since we will take the renormalization scale to 
be Qo = 500 GeV, the running of the gauge couplings from Mz to this scale is negligible. 
We also omit the couplings g\ and g 2 in the calculation of the one-loop effective potential. 
However, we retain the terms with these couplings at tree-level since they include the quartic 
Higgs scalar interactions responsible for the non-zero Higgs vev of the Standard Model (SM). 



III. THE EFFECTIVE POTENTIAL 



Quantum corrections may affect whether spontaneous symmetry breaking occurs in a 
field theory. The vevs of the scalar fields are the values of the classical fields at the minimum 
of the effective potential. We use the one-loop correction to the effective potential in the 



dimensional reduction, DR, renormalization scheme pA 

1 
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(3) 



where the supertrace, Str, is over all color and spin degrees of freedom with a minus sign 
for fermions. Q 2 is the square of the renormalization scale, which we take to be equal to 
Qq. Ai 2 is the field dependent mass matrix; for instance, for a theory with scalar fields 
6i,i — 1, . . . ,n 



Q2y 

Mi 



(4) 



M. as well as V e ff are function of the classical fields (<j>i) c = (O|0«|O). Since the top Yukawa 
coupling, ht, and the strong gauge coupling, #3, are large one may expect significant contri- 
butions to the effective potential from one-loop corrections. The one-loop effective potential 



is also more stable under a change of renormalization scale [pl}| . Since the parameters ap- 
pearing in the equation for the one-loop correction to V e ff are the renormalized parameters 
we consider the resulting CCB constraints as being limits on these renormalized parameters 



in the MS scheme at momentum scale Qq. If large logs of ratios of massive parameters 
appear in the effective potential, such as In (jn 2 /m 2 ^j or In (d) 2 /m^j, then one may, in prin- 
ciple, sum these using the renormalization group equation |T5 ]. However, this is difficult 
when there are multiple scalar fields and masses. We only consider masses, m^, and the 
renormalization scale, Qq, which differ by less than two orders of magnitude. Furthermore, 
since we are only interested in the effective potential near its minimum where (<f>i) c ~ Trij, 
i.e., the classical field near the minimum of V e ff is of the same order of magnitude as the 
masses, the logarithms appearing in V e ff are small. Renormalization group improvement of 
the effective potential in this case is unnecessary. 

We also assume that the Lagrangian parameters are real. In a field basis where the Higgs 
fields are real and positive and ht and the scalar masses, mf, real, only one of the parameters 
At, fi, or m\ may be made real by redefining the fields. However, it was shown in Ref. (16| 
that the presence of complex phases in these parameters, greater than about 10 -2 , gives too 
large of a contribution to the neutron's electric dipole moment due to gluino loops. 
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IV. NUMERICAL CALCULATION 



We evaluate the one-loop effective potential using Eq. (]5p by diagonalizing the mass 
matrix M.. If one chooses the Landau gauge for the gluons, in which the propagator is 



— i 



Qpiv k^k y 

k 2 



k 2 + le 



(5) 



the generalized mass matrix Ai is block diagonal in the gluon fields, bosonic fields, and 
fermionic fields. The top squark fields may also be rotated by a global SU (3) transformation 
so that they are of the form 



/ \ 

o 

V J 



tlR + lt\I 



V 



t2R 




(6) 



/ 



with qui, tm, tu, and t 2 R are real. We then use only the real part of the one-loop effective 
potential; the imaginary part is related to the decay rate of the vacuum as shown in Ref. 

The global minimum of the effective potential is found by calculating the local minimum 
using a standard algorithm from Ref. |18|] starting from a set of field values on a rectangular 
grid. If this global minimum occurs for non-zero values of the squark fields then there is 
a CCB vacuum. This process is repeated for a set of soft-breaking parameter values on a 
rectangular grid in parameter space. A quadratic surface is then fit to the parameter points 
on the boundary between those that give a symmetric vacuum and those that lead to a CCB 
vacuum. The surface is fit to the boundary points by varying the coefficients in the equation 
of the surface, (Pi, . . . ,P%), until the average of the distance squared from the boundary 
points to the nearest point on the surface, (d 2 ), is minimized. The equation of the surface is 



A\ + P X A\ + P 2 m\ + P 3 m 2 2 + P 4 mi + P h m\ + P^n\ + P 7 m\ - P 8 = 



(7) 



with (At, m\, m 2 -, m~) a point in parameter space. These soft-breaking parameters are defined 
in the appendix. This particular parameterization of the surface is chosen solely because 



of its simplicity, but as will be seen later it gives an accurate characterization of the CCB 
boundary. 

Because of the large dimension of the parameter space in which we want to find CCB 
bounds, it is necessary to initially constrain some of the parameters values in order to reduce 
the calculation time. Two of the soft-breaking parameters, m\ and m 2 , are constrained using 
two different methods. In method S the values of these parameters are chosen such that 
there is a Standard Model minimum in the effective potential with (q) = (t c ) = and 
(Hi) = v cos/3, (H2) = v sin/3 with v = 246 GeV. This is done by solving the analytical 
expression for the one-loop effective potential with zero vevs for the top squarks in terms of 
these parameters. These relations are shown in Eq. fl2"2"|). In the other procedure, method 
F, we simply fix m\ and m\ at constant values, i.e., independent of the other soft-breaking 
parameters. 



V. RESULTS 



A. CCB Bounds for a Simplified Model 



We begin with a numerical analysis of the model with three real scalar fields H, Q, and 
T of Ref . IJ . The potential in this model is 



h 2 V = (Q 2 T 2 + Q 2 H 2 + T 2 H 2 ) - 2AQTH + m 2 Q Q 2 (8) 

1 

8h 2 



+ m 2 T T 2 + m 2 H H 2 + ^ 



gl{H 2 + l -Q 2 - A -T 2 ) 2 



+ g 2 2 (H 2 -Q 2 ) 2 + ^g 2 3 (Q 2 -T 2 ) 2 

This is a simplified version of the potential described in the appendix with Q = q, T = t c , 
and H = H°. The CCB bound 

A 2 < 3(ml + m 2 Q + m^), (9) 

derived in Ref. []5| follows from minimizing the potential only in the equal-field direction, 
i.e., Q = T = H. This is valid, in general, only if the D-term dominates, or if g%jh ^> 1 



with gi the smallest gauge coupling. 

We perform a numerical analysis of the potential in Eq. (|8|) using the parameter value 
ranges shown in Table | and with h = 0.1. The gauge couplings are fixed at their running 
values at the weak scale Mz as stated earlier. The resulting best fit to the CCB boundary 
is shown in Table |TJ. Clearly, this is close to the analytical result of Eq. as it should be 
with a small value for h. 

Next we repeat the analysis for the potential of Eq. (H) with the same soft-breaking 
parameter ranges of Table |, but we set h = 1.0. This is the proper value for the top quark 
Yukawa coupling for small tan/?. The corresponding CCB bound is given in Table |I[ 

The CCB bound for h = 1.0 is more stringent than the bound for h = 0.1, i.e., for a given 
set of parameter values within the range given in Table |, (mf, m|, ra~), the value of the A 
parameter at which the vacuum becomes CCB is lower. The fractional difference in the A 
parameter bound between these cases is largest when the other parameters, (m|,m~,m~), 
are near the lower end of their range in Table Q, where A A/ A « 15. Also the value of A 
at the CCB bound in this case is small. The average fractional difference over the entire 
parameter ranges shown in Table | for the two bounds is AA/A m 0.14 

B. CCB Bound for the MSSM 

Next we examine the CCB bounds obtained from of our approximation to the MSSM 
with the Lagrangian given in the appendix. In all of the following analyses the mass units 
are TeV. We also fix tan/? = 1.5 and fi = 0.4. The renormalization scale for the one-loop 
calculation of the effective potential is set at Q — 0.5. 

We first examine the bound for method F with m\ = 0.65 and m\ = 0.40. These masses 
result in a value of m\ = 0.26 if a S.M. minimum were required, i.e., zero squark vevs. 
These values for m\ and m\ were chosen since they are consistent with a S.M. minimum 
for tan/5 = 1.5 and they are close to the renormalization scale, Q = 0.5. The ranges of 
parameter values for which the effective potential was calculated are shown in Table [JTI]. 



S 



The CCB bounds for the analyses of both the tree-level and the one-loop effective potential 
are shown in Table |V|. Parameter values which give a potential that is unbounded below 
are not used in determining the best fit surface for the CCB bound. This includes values 
with m 2 < 0.15 for the tree-level potential since in this case the potential is unbounded for 
m\ < 2m\ — m\ and zero squark vevs. We do not examine negative m 2 values since in this 
case, for the small range of m~ and m~ values for which the potential is bounded below, 
the vevs for the CCB. vacuum are too small for the numerical methods to distinguish this 
minimum from the minimum with all vevs vanishing. However, according to Eq. (|2ll) , the 
Higgs soft-breaking mass, m 2 H2 , is negative over part of the parameter range examined. 
Next we calculate the CCB bounds using method S, in which m 2 and m\ are fixed by 



requiring a S.M. minimum. The parameter values listed in Table [III] were used in determining 
the bounds. Since there is no Standard Model minimum for m 2 , < (60 GeV) 2 , we do not 
examine m 2 , values in this range. The CCB bounds calculated by finding the global minimum 
of the both the tree-level potential and the one-loop effective potential are given in Table 



VI. DISCUSSION 

The quadratic surface of Eq. ([?]) is sufficient to provide an accurate characterization of 
the numerical CCB bound. All of the least-square fits give the average distance squared, 
(d 2 ) < 1.0 x 10 -2 . The one-loop CCB bound calculations yield a larger value for (d 2 ) because 
the longer calculation time requires using less parameter grid points than for the tree-level 
calculation. Typically it takes around 20 days of cpu time on a Digital Alpha workstation 
to calculate the one-loop CCB boundary points and perform a least-squares fit so we are 
limited by the computation time. 

As stated in section |V A] , the numerical CCB bound for the simplified Lagrangian of 
Ref. U with Yukawa coupling h — 1.0 is significantly different from the numerical bound 
with h = 0.1, and hence also different from the analytical bound of Eq. when the other 
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soft-breaking parameters are small. However, for the remainder of the parameter values 
tested, the numerical bound with h — 1.0 is quite close to the analytical bound. One would 
not expect the CCB bound of Eq. to be correct for large h. One possible explanation is 
that the g\ D-term is large enough at least to insure that the minimum is in the direction 
Q = T. Also for h — 1.0 the numerical bound does give a more stringent CCB bound than 
that of Eq. (||) over the entire range of parameters tested. 

We present some contour plots of the CCB bounds with m\ = 0.25 for comparison. 
The contours show the value of A t on the CCB boundary, i.e., lower values of A t result in 
a symmetric vacuum whereas higher values result in a CCB vacuum with nonzero squark 
vevs. Figs, p] and ^| show the bounds for method F using the tree-level and one-loop effective 
potential respectively. Figs. ^| and f| show the tree- level and one-loop method S CCB bounds. 
Fig. |5| shows the analytical bound of Eq. (|D with m\ = 0.25. 

Since there is not the additional constraint of requiring a Standard Model minimum for 
the Higgs field for both the method F CCB bounds and for the analytical bound, one may 
compare these CCB bounds. The tree-level CCB bound on the A t parameter for method 
F is lower than the analytical bound of Eq. (|9|) for the entire range of parameter values, 
(m?,, m^,mf), shown in Table fTJ. The one-loop correction to the effective potential raises 
the value of A t for the CCB bound over about 70% of the parameter range considered. 
However, even with the one-loop corrections, the CCB bound for the MSSM potential is 
more stringent, i.e., gives a lower value for A t , than the bound of Eq. (||) for > 95% of the 
parameter range. 

The one-loop corrections for the CCB bound calculated using method S give a lower A t 
value than the tree-level bound over the entire range of parameters examined. For values of 
the parameters, (m^, m|, that give a small value for the A t parameter CCB bound, the 
one-loop corrections to the effective potential makes the CCB bound significantly stricter. 

In conclusion, CCB bounds on the soft-breaking parameters of the Higgs and top 
quark/squark sectors of the MSSM provide important constraints for these parameters. 
These constraints may be expressed as a maximum value of the A t parameter for given 
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values of the remaining soft-breaking parameters. The numerical CCB constraints that we 
calculated give more stringent CCB bounds than the analytical constraint of Eq. (|9]) for 
most of the ranges of parameter values considered. Because of the large top Yukawa cou- 
pling, one-loop corrections to the effective potential may result in significantly different CCB 
bounds than those for the tree-level potential. 
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APPENDIX 

Tree-level MSSM Potential 

The potential may be divided into several parts. A sum over group and spinor indices, 
where applicable, is implied. The supersymmetric D-terms are 

Vd = \gl {-HlH, + H\H 2 + i|g| 2 - t\*f) (10) 
1 3 

+ q9% E + H\r a H 2 + q*r a q 



u a=l 
1 a=l 

where gi, g 2 , and g 3 are respectively the U(1)y, SU(2), and SU(3) couplings, r a are the 
Pauli matrices and T a are the antihermitian generators of £77(3). Using the relations 

T^flx = 28uSjk - 5ij5 k i (11) 

and 



T ij T ki = j (tiatjk ~ ^ijtki ) (12) 
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the SU (2) contribution becomes 

1 

VSU(2) = g#2 



HlH,) 2 + (H\H 2 f + (gtg)^ 
- 2 ([H\Hi)(H\H 2 ) + (HlHjtfq) + (H}H 2 )tfq 
+ A(\H\H 2 \ 2 + \H 1 q\ 2 + \Hlq\ 2 ) 



(13) 



and the SU (3) one is 



V S U(3) = 



l^q) 2 + ^~ c ) 2 - (#) f (#) + ^SX*^) 



(14) 



The superpotential or F term is 



V F = h 2 (\q\ 2 \H° 2 \ 2 + \tf\H° 2 \ 2 + |#| 2 ) + htuffiH?) + /i.e., 



(15) 



with h.c. denoting the Hermitian conjugate and H® and H 2 are the neutral components of 
the Higgs scalar doublets. The Higgs scalar and fermion doublet components are 



Ho 



(16) 



The quark-squark-gluino interaction terms are 



V q g X = tV2g 3 (tP L X {a) T a t c * - t c T a X ia) P R t + q*T a XP L t - tT a qP R X {a) ) 



(17) 



Pl,r are the projection operators for left- and right-handed chiral spinors, |(1 ±7 5 ), t is the 
four component spinor field for the top quark, and A*") are the Majorana gluino fields. The 
quark- squark-Higgsino terms are 



V q . 5 = //, ( t°H 2 P L t + t?*tP R H° - qll',11':! - q IIJ> 



(18) 



The Higgsino interaction terms are 



V m = fjL (H°H° - H~H+) + h.c. 



(19) 



Finally the SUSY soft-breaking terms are 
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V soft-breaking — TH^H^i + W? H2 H\H2 — ^3^06-^1-^2 

+ m?|g| 2 + mfltf + im A A (a) A (a) + h t A t qt c H° 2 + h.c. 



(20) 



With the addition of the supersymmetric Higgs interactions, the masses for H\ and H2 
become 



m H1 + fx 



(21) 



m 2 = m^ 2 + fi 



respectively. 



One-loop Effective Potential with Zero Squark VEVs 



If the squark vevs are zero the one-loop contributions to the Higgs effective potential 



from top squark and quark loops may be written in an analytical form jnj. After requiring 
that the minimum of the effective potential be at H\ = v\ and H2 = v 2 and solving for m\ 



and m§ we obtain 



m 2 = (777.3 — m l cos 2/3) tan (3 



{ht A t (A t tan f3 + fi) 



(22) 



16vr 2 



f(m'i) - /(mi 2 



m 9 



mi 



+ ^ 2 tan/3(/(mi 2 ) + /(m 2 2 ) - 2/(m 2 )) 
m 2 = m 2 tan (3 — m 2 z cos 2(3 



+ 



167T 



■h\ii(A t tan (3 + li) 



f{jn{) - f(m 2 2 ) 

"1 "2 
m 2 — mi 



with f(m 2 ) = m 2 (ln(m 2 /Q 2 ) — 1). The definition tan (3 = t^/Vi and the tree-level Z boson 
mass, m% = {vf + v^igf + g^)/^, were used. The tree-level relation follows by including only 



the first term in the above equations for m\ and m\ 
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TABLES 

TABLE I. Grid of parameter values used to calculate the CCB bound for the Lagrangian of 
Eq. d). 



Parameter Min. Value Max. Value Grid Spacing 

A 0.0 3.0 0.1 

m 2 2 -0.250 0.750 0.075 

m 2 Q 0.0 0.750 0.075 

m\ 0.0 0.750 0.075 



TABLE II. CCB bounds for the potential of Eq. (||) with h = 0.1 and h = 1.0. The coefficients, 
Pi, are defined in Eq. ([?]). 



Parameter 


h = 0.1 


h = 1.0 


Pi 


0.7246 


2.723 


P 2 


0.1608 


3.527 




-3.808 


-8.310 


Pa 


0.1991 


3.370 


P 5 


-3.830 


-8.228 


Pe 


0.1985 


3.452 


Pi 


-3.836 


-8.310 


Ps 


0.4911 


0.5877 


(d 2 ) 


2.8 x 10" 4 


4.1 x 10 
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TABLE III. Grid of parameter values used to calculate the tree-level and one-loop CCB bounds 
for the MSSM Lagrangian given in the appendix. 



Parameter Min. Value Max. Value Tree-level Grid Spacing One-loop Grid Spacing 
A t 0.0 3.0 0.1 0.4 
m\ 2.5 x 10~ 3 0.7525 0.125 0.25 
mf 2.5 x 10~ 3 0.7525 0.125 0.25 
mf 2.5 x 10" 3 0.7525 0.125 0.25 



TABLE IV. Tree-level and one-loop CCB bounds for the MSSM Lagrangian using method S. 
The coefficients, Pi, are defined in Eq. (0). 



Parameter 


tree-level 


one-loop 


Pi 


6.192 x 10 3 


5.370 x 10 1 


P 2 


3.078 x 10 3 


1.203 x 10 1 


Pz 


-3.107 x 10 3 


-4.913 x 10° 


Pa 


6.721 x 10 3 


6.872 x 10 1 


P 5 


-1.162 x 10 4 


-1.120 x 10 2 


Pe 


6.857 x 10 3 


7.268 x 10 1 


Pi 


-1.175 x 10 4 


-1.068 x 10 2 


Ps 


-1.305 x 10 3 


-1.233 x 10 1 


(d 2 ) 


1.2 x 10~ 3 


4.8 x 10~ 3 
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TABLE V. Tree-level and one-loop CCB bounds for the MSSM Lagrangian using method F. 
The coefficients, Pi, are defined in Eq. (^). 



Parameter 



tree-level 



one-loop 



Pi 


4.431 x 10 2 


3.000 x 10° 


p 2 


9.397 x 10 2 


-1.307 x 10 1 


P3 


-1.452 x 10 3 


5.948 x 10° 


Pa 


3.478 x 10 2 


7.767 x 10° 


P 5 


-6.650 x 10 2 


-9.086 x 10° 


Pe 


3.553 x 10 2 


-7.102 x 10° 


Pi 


-6.700 x 10 2 


3.374 x 10° 


Ps 


-2.346 x 10 2 


-2.833 x 10" 


(d 2 ) 


1.6 x 10~ 3 


9.4 x 10~ 3 
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FIG. 1. Tree-level CCB bound for method F given in Table fv|. m| = 0.25 and each contour is 
labeled with the corresponding maximum A value. 
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FIG. 2. One-loop CCB bound for method F given in Table [V|. 
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FIG. 3. Tree-level CCB bound for method S given in Table IV. m| = 0.25 and each contour is 
labeled with the corresponding maximum A value. 
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FIG. 4. One-loop CCB bound for method S given in Table IV. m| = 0.25 and each contour is 
labeled with the corresponding maximum A value. 
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FIG. 5. Analytical CCB bound for the potential of Eq. (g) given in Eq. @. = 0.25 and 
each contour is labeled with the corresponding A value. 
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